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Abstract In this paper we refer to the reconstruction formulas given in An- 
dersson's On the determination of a function from spherical averages, which 
are often used in applications such as SAR 1 and SONAR 2 . We demonstrate 
that the first one of these formulas does not converge given physically rea- 
sonable assumptions. An alternative is proposed and it is shown that the 
second reconstruction formula is well-defined but might be difficult to com- 
pute numerically. 

1 Introduction 

The determination of a function from spherical averages is a problem often 
encountered in physical applications such as SAR and SONAR. The work 
related to this topic, which has lead to a great amount of insight and refinement 
today, began with the proposal of a reconstruction formula by Fawcett 0. The 
mathematical analysis of the problem was later improved [2| and two refined 
reconstruction formulas were derived. This sparked a host of activity that 
branched to several areas 0, 0, 0, 0, so that most of today's research is 
based on Anders sons 's ideas. 

In Q two reconstruction formulas were derived from the fourier inversion for- 
mula. But it was neglected to examine whether they are properly defined. In the 
following it will be shown that one of these is useless, if the function to be re- 
constructed has physically sensible properties, because of a divergent integral in 
this case. An alternative will be presented and it will be shown that the the other 
reconstruction formula is well defined but might me difficult to compute numeri- 
cally. 

At first, for the benefit of the reader, some results from will be recalled. Note 
that besides the aforementioned problems, a few minor errors occurred, which do 



1 SAR: Synthetic Aperture Radar 
2 SONAR: SOund NAvigation and Ranging 
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not essentially obscure the results in [2|. For a detailed analysis of these errors see 

ED- 

1.1 Definition 

Let fiGN. Then 

1. y(R n+1 ) is the Schwartz Space. 

2. y e (R n+1 ) := {ip E y(R n+1 ) : (p{x,-y) = ip{x,y),Vx E R n , y E R}. 

3. ^ r (R n x R n+1 ) := {</? G J^(R 2rt+1 )|VorthonormaJ transformations 

U : R n+1 -> R n+1 Vx E R n , z E R n+1 : <p(x, z) = (p(x, Uz)}. 

4. y^(R n+l ) and J^f(R n x R n+1 ) are the dual spaces of y e (R n+l ) and 
y r {R n x R n+1 ) respectively. 

5. Let f E y e (R n+1 ). Then the operator R is defined by 

Rf(x,r) := J f(x + r£,rr))dS n (£,r]). 
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6. S n denotes the unit sphere in R 

For simplicity in the following sometimes only y, y e , etc. is written instead of 

y{R n+l ), y e (R n+l ), etc. 

1.2 Remark 

It is easily seen that f E y e and g E y r , if f E y e and g E y r , respectively. 
The essential result in Q is the (Fourier-) inversion formula: 

1.3 Theorem 

Ify e (R n+1 ) is given the topology ofy'(R n+1 ), the mapping 

R : y e {R n+1 ) -> y;(R n x R n+1 ) 
is continuous and can, by continuity, be extended to a mapping 

r : y^{R n+1 ) -> y;(R n x R n+1 ). 

The range of this extended mapping R is the closed subspace 

y'(R n x R" +1 ) = {QE y' r {R n x R" +1 ) : supp g C {(^r,) : \\ v \\ > ||£||}}. 



R is one-to-one and the inverse mapping 



R- 1 : ^ cone (R" x - ^(R n+l ) 



is continuous. Moreover, if g = Rf and iff(£, rj) or g(£, rj) are integrable for £, r] 
in some open set, then 




or 



1 1ST- 



v\(M\\ 2 + v 2 )^m vfF+?)- 



(2ir) n 2 



respectively. 



Proof: 

El Theorem 2.1] 



□ 



In this theorem the operator R describes in the case of SAR and SONAR the 
measurement of the reflectivity function / that represents the ground reflectivity. 
The measurement is modelled as a 5-impulse wave that propagates as concentric 
spheres. The ground is approximated as a plane. The single scatter approximation 
for the wave hitting the ground results in integrals over circles. The Fourier trans- 
form g of the data g can under certain conditions be used to extract the Fourier- 
transform / of the reflectivity function /. 

1.4 Remark 

Let f G Y e (R n+1 ). Then for g = Rf e ^'(IET x R n+1 ), x G W 1 and r > in 
the following g(x, r) is sometimes written with abuse of notation. This is justified 
because g depends only radially on the last n + 1 variables. 



In Q two reformulations of this formula were proposed, of which the first was 
already essentially given in CJ. 

1.6 Corollary 



1.5 Definition 

Let g G y r {R n x R n+1 ). Then 




With c. 



1 \s' 



two reformulations of the inversion formula are possible: 



(271-)™ 2 
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1. Forg G y r 

d 



f = c n Hy—A~R*g 



'dy 

with the Hilbert transform in y, H y , and the Laplace-Operator A = A x + 
■ This formula is essentially also given by Fawcett OF- 

2. Forg G y r , cone : = {g G Y r : supp g C {(£,77) : |M| > ||£||}} 

f = c n R*Kg 

with the pseudodifferential operator K defined hy Kg(^,rj) = 
Proof: 

El Section 3] 

□ 

The corollary states that under certain restrictions it is possible to reconstruct the 
reflectivity function / directly from the data g without taking the detour through 
the Fourier space. 

1.7 Remark 

Note that the essential restriction of this formulation of the corollary in compari- 
son to 0/ is that the data g has to be in y r or S r>cone respectively. This is necessary, 
because otherwise the application of R* to g in the first case is undefined or the 
application of R* to Kg in the second case. 

Unfortunately g = Rf is usually not in 5^ T . Therefore the two reformulations of 
the inversion formula are only valid in the distributional sense with an appropri- 
ately defined R*. 

For example a physically reasonable f G y e , f : R 2 — > [0, 00) with f(x , y ) > 
for some x , y G R yields f(x, y) > c > for all (x, y) G K e (x , y ) with ap- 
propriate c, e > 0. 
Moreover 

2ng(z, y/(z -x) 2 + y 2 ) = \ =^ / f(z + r 1 ,r 2 )da(r) 

yj{z — xy + y z J 

ll r ll=\/ (z—x) 2 +y 2 

with r = (n, r 2 ). For \/ [z — x ) 2 + y% > e we obtain with a simple geometrical 
consideration: 



2ng(z, yj(z- x ) 2 + y 2 ) > 



ec 



2^{z-x Q ) 2 + y 2 
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2nR* 



*g(x ,yo) = / g ( z, yj \z - x ) 2 + y%\ dz 



TR 



ec f 1 

> — / — ; dz = OO. 

~ 2 J yfc ~ X ) 2 + Vl 

Therefore the reconstruction of a non-negative function f with f(xo, yo) > for 
some (xq, yo) e R 2 is impossible using the first reconstruction formula. 

This result is in accordance with a result from Nessibi, Rachdi and Trimeche [9]. 
They gave reconstruction formulas for functions / with 

oo 

J P(r)f(r, x) dr = 
o 

for all x G R™ and for all one-variable polynomials P. 



2 Properties of the function g = Rf 

Before showing an important property of the data function g = Rf, which will be 
necessary for the derivation of the new reconstruction formula, some definitions 
and a result from [2| will be needed. 

2.1 Definition 

Let fiGN. Then 

1. ^(R n+1 ) := C^°(R n+1 ). 

2. := {ip e @(R n+1 ) : cp(x, -y) = ip(x,y),Vx ER n ,yE R}. 

2.2 Corollary 

Iff E y e (R n+1 ), then g = Rf G ^ r '(R n x R n+1 ). 
Proof: 

El Section 2] 

□ 

2.3 Proposition 

Iff G y e (R n+1 ), then g = Rf G 

Proof: 

The change of differentiation and integration is justified, because / G ^ e (R n+1 ). 

□ 
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2.4 Theorem 

Iff E y e {n n+1 ), then g = Rfe L X (R" x 



Proof: 

With fey, also / eyci 1 . Theorem Q1 implies that 







M NrviMl 2 -ll£ll : 



otherwise. 



1^,^)1^= (2tt)"— - 
The substitution p' = \\r)\\ results in 



l>ll£ll 



m,y/\\v\\ 2 -U¥ 



IN n -VIN 2 -lieiP 



d^drj. 




= 2(27r) n 

E" p'> 

The substitution p = p' + ||£|| leads to 



p'm, Vp 12 - lien 2 ) 



Vp' 2 -CT 



d^dp' . 



w// 



R" p>0 



(p+ll^ll)/(e,v / P 2 + 2p||ell) 



/ is in =y, therefore 

< 2(2vr) n 




VP 2 + 2p||ell 



C(P +||eil) 



d^dp. 



R™ p>0 

= 2(2vr) n 
< 2(27r) n 



c(p+neii) 



d£dp 




R n p>0 




R n p>0 



v /p2 + 2 p||e||(l + ||e||+p)" +3 

C 

Vp 2 + 2p||e||(l + ||^||+p)" +2 



d^dp 
d£dp 



and the integral 




C 



R" p>l 

converges. So only the integral 



^ p 2 + 2p u\\(i + U\\ + p) n+2 



d£dp 




C 



R n 0<p<l 



Vp 2 + 2p||e||(l + ||ell+p) 



n+2 



d£dp 



remains to be examined. 

C 



IS 



— dfdp 

+ IKII +/>)"" 

-/„ X TSoTra^^ 



The integral 



^W(l+llelir 2 



V v^(i + lieil) n+2 ' 



>i 



converges. For the proof it therefore suffices to show the existence of 



4<i v^( i+ n«ii) n+2 



di 



+2 



IKi 



v^(i + lieil) n 
f 2C dr 

ll€ll<i 

The substitution r = ||£|| yields 

= 2C\S n - 1 \ I r n ~$dr < oo, 

because n > 1. 
2.5 Remark 

An analogous proof shows rjg G L 1 . 



□ 



3 Modifications to Andersson's first inversion for- 
mula 

3.1 Definition and properties of a modified R* 

Now a modified version of the operator R* is introduced and its properties are 
discussed. 
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3.1 Definition 

For f e 3> e andg = Rf e C°°(R n x K™ +1 ) we define 




x — z\\ 2 + y 2 ) dz. 



This slight modification by an additional derivation turns out to ensure the conver- 
gence of the integral applied by the operator R* d under the minor and physically 
feasible constraint that / is in 3> e . Thereby the formulation of a mathematically 
exact reconstruction formula is possible. 

3.2 Corollary 

R* d g is well-defined. 

Proof: 

Let x G R™, 1/6R. Then 



/ G 9i e and therefore it is possible to exchange differentiation and integration. 




x - z\\ 2 +y 2 ,r]^\ 



x - z\\ 2 + y 2 )dS n (^r])dz . 



x - z\\ 2 + y 2 ,vV\\ 



x - z\\ 2 + y 2 ) dS n (£,,r])dz 




= Q-(v/)(, + Mf 



x - z\\ 2 + y 2 ,r]>/\\ 



x — z 



2 + y 2 ) dS n (£,r))dz 




fll(V/)(«+eVP - ^ll 2 + i^WP - ^ll 2 + 2/ 2 )ll dS n {Z,rj)dz. 



With r = (£, 77) a/H^ — z\\ 2 + this is 




\\(Vf)((z,0)+r)\\d<j(r)dz. 
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fe@=>\(R* a g)(x,y)\ 



\y\ 



(^/\\x- z\\ 2 + y 2 ) n+1 



max (|| V/ll) diam(supp /) dz < oo 



□ 

3.3 Corollary 

For f G @ e and g = Rf, R* a g G . 
Proof: 

This is guaranteed by the estimate in Corollary 13 .21 

□ 

3.4 Definition 

Letfe y'(R n ) and if e y(K n ). Then 

< f,<P >j?(b»)'= 

is the functional f applied to the test function ip. Here the subscript y{K n ) is a 
reminder of the space of the test function (p. 

Now analogously to [2] an expression for R* d g is derived. 

3.5 Theorem 

Let fe3> e andg = Rf. Then R*g(£, rj) = irjg(^ >/||£|| 2 + r? 2 ). 
Proof: 

Let ip G y, f G % and g = Rf. Then 



< R%9, ¥> >^(E.»xB.)=< Rig, >y(R» X R)= J (Rd9)(x',y)<p(x',y) dx'dy 

E"xR 

= J J J-S'fo V\\ x> - z \\ 2 + V 2 ) dz<p(x', y) dx'dy. 
As R* d g G 5?' and <p G y, Fubini's theorem implies 

= j j J-tffo y/\\x' - z\\ 2 + y 2 )<p(x', y) dx'dydz. 

R™ R"xR 

Substituting x = x' + 2; we obtain 

= j J -jf9(z, V\\ x \\ 2 + y 2 )<p(x + z, y) dxdydz. 

K.™ R"xR 



9 



Letijj(x,y, z) := ip(x + z,y). Then ip(x, y, z) e y(R n xRx R") and this leads 
to 

Q 

£fy~9(z, V\\ x \\ 2 + y 2 )i ; ( x , V, z ) dxdydz 

R n E"xE 

d . IT ~ . ~ 

=< W >y(R"xRxE»)=< dy 3 >y(rxRxR»)=< ^g-iW >y(R"xRxR») 



ft" 

iV9((,VM\\ 2 + V 2 ) J J (27r)- 2 "- 1 e i «^«> + ^ +<z ^»^(a;, 2 /,z)rfa:rf^C^ 

'xIR" R"R"xR 

j J ir]g((, VU\\ 2 + V 2 ) J J (2n)~ 2n - 1 e t ( <x 'S >+ y r i +<z '< > ^(x+z,y)dxdydzd(d£dri. 



R"xRR" 



R"xlE" R"R"xR 

As rjg G L 1 and tp e S fi , Fubini's theorem gives 



iri9((, VUW 2 + ^ 2 )(27r)- 2 "- 1 e l(<:E ' ?>+OT+<z ^ > V(a:+2, y) d(d£dridxdydz. 



R" R n X R R n X R R™ 

The substitution x' = x — z yields 



irjg 



(C, Vll^ll 2 + ^ 2 )(27r)- 2 "- 1 e l(<:E '' f>+OT+<z ' c - 5> V(^, y) dCdidrfdx'dydz 



R™ R n X R R n X R R™ 



R n R™ x R R n 

Continuing in the distributional sense we obtain 



*vg(C, vfTF+?)^(( - £M£, v) d(dtd v 



R n x R R n 



ivgit, y/U\\ 2 + V 2 MZ,v)dtdr, 

R"xR 

=< ivgi, V >y(R"xR) • 
Here g\(£, 77) := <?(£, ^U\\ 2 + v 2 ) and with this definition g\ e J^'(R n x R). 

□ 
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3.2 A modified inversion formula 

With the results of the preceding sections a well-defined reconstruction formula is 
attainable also for physically meaningful reflectivity functions. 

3.6 Theorem 

71 — 1 

Let f G 3) e and g = Rf. Then f = c n H y A^~ R* d g with the constant c n := 

1 I S n I 

j^y; l y 1 and the Hilbert-transform H y . 
Proof: 

Let / e @ e . It follows from Theorem 12 .4! and Theorem II .3! that 

f = f = CndvKUW 2 + v 2 )^ g)~ 

= c n {-i sgn(r])(\\£\\ 2 + n 2 )^ ingj. 

Theorem 13.51 yields 

= c n {-isgn{r 1 )m\ 2 + V 2 )^R^9y 
= c n {-i sgn(r])(A^R* d g)J 
= c n ((H y A^ Rlg)J 

71 1 

= c n H y A — R* d g. 

□ 



4 Well-definedness of Andersson's second inversion 
formula 

In the following it is shown that the second reformulation of the inversion formula 
in Q is well-defined. To this end at first some properties of Kg are derived. 

4.1 Proposition 

Let f ey e andg = Rf. Then Kg e L p (R n x R n+1 ) for all p > 0. In general g 
is only in C°° and not for example in L 2 . Therefore g has to be computed in the 
distributional sense. 

Proof: 

Theorem 1 1 .31 implies that 

9(Z,v) = < \s n \ iMi»-yi|7/ll 2 - neil 2 

I otherwise 
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2 



(^/(f. VIRFIF) for |M| > Ml 
otherwise 



2 



<<; ( 27r ) n ^( 1 + v / iieii 2 + ihii 2 -neii 2 i)-^ forhn>neii, 

otherwise 

2 2n + 2 

<<; ( 2 *) n j^( 1 + Nir~ for|M|>||a 

otherwise 

2 2n + 2 

< (270^(1 + max(||aiM|)r—. 



4.2 Corollary 

Let / G ^ e with /(0, 0) ^ and # = Rf. Then g,g£ L 2 (R n x R n+1 ). 
Proof: 

Theorem 11.31 yields 



□ 



(2l) ^/yra forwi>iia 

otherwise 

|s(£,J7)| d£dr) = {2n) —— / , „ 2 „ ,,,, ,| 2 rr-^d^dr). 



\Sn\2 J 11^11^-2(11^12 
R"xR»+ 1 IMI>II£II 

The substitution ||n|| = p' results in 



^ W\l J ?- M' ^ 

H» p>>U\\ 



The substitution p' = p + ||£|| leads to 



(2 „ 4 /• /■ (p+||ell) 2 -1/(e, v / P 2 + 2p||ell)| 2 



1 ; i^i y y P 2 +2piKn s ^ 



E" p>0 



We consider 

(p+II£II) 2 -1/(£, Vp 2 + 2pU\\)\ 2 



P 2 + 2pU\\ 

<2°<P<2 

12 



d^dp 



and assume without loss of generality /(£, y/ p 2 + 2p||£||) > 1 for ||£|| < |, < 
P< l r 

- y y ttw^' 

uen<|o<p<| 

With the substitution r = ||£|| we obtain 

0<r<i 0<p<^ 



For n = 1 this is 



> / / ^±I- d rdp> J J ±drdp = oc. 



0<r<±0<p<± 0<r<iO<p<i 

For n > 2 this results in 



> / — drdp = oo. 

-J J p 2 + 2pr H 

0<r<| 0<p<\ 

□ 

4.3 Remark 

Unfortunately /(0, 0) 7^ for a reflectivity function f with physically feasible 
properties. Moreover g is not continuous and therefore g L 1 , so it might he 
difficult to compute g numerically with sufficient accuracy. 

4.4 Theorem 

Let / 6 y e and g = Rf. Then R*Kg is well-defined with g computed in the 
distributional sense. Further 

Proof: 

It follows from Proposition ^, ll that Kg £ L 1 fl L 2 . Therefore 



Kg(Z, VU\\ 2 + v 2 ) 



e 



R" 



A< ^ Z> J e-^'^+^Kgiz, y/\\x\\ 2 + y 2 ) dxdy J dz 

R"xR / 
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: J I e -i<^> J e - l ^' x '~ z>+vy) Kg(z, y/\\x'-z\\ 2 + y 2 ) dx'dy J dz 



e -i<Z,z> e -H<eJ-z>-hw) Kg ( Zi v /||. x '_^||2 + y 2) dzdx ' dy 



R"xEE" 



R"xE R n 

v s- 



=R*Kg(x>,y) 

= R*K~g(£,ri) 

and the change in the order of integration and the existence of R*Kg is justified 
by Fubini's Theorem. 

□ 

With this result it can be shown that the second reformulation of the inversion 
formula in Q is well-defined. 

4.5 Corollary 

Let /ey e and g = Rf. Then f(x, y) = c n R*Kg(x, y) with a constant c n . 
Proof: 

Theorem 1 1 .31 implies that 

m,v) = CnlvKUW 2 + r?)^g{t, vWT?) 
with a constant c n . 

=► m Tj) = c n Kg{i, 

= cjFKg{^rj). 
This completes the proof, because / is in S e . 

□ 
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